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ABSTRACT 

Recent theoretical work has suggested that Ly« nebulae could be substan- 
tially polarized in the Lya emission line, depending on the geometry, kinematics, 
and powering mechanism at work. Polarization observations can therefore pro- 
vide a useful constraint on the source of ionization in these systems. In this 
Letter, we present the first Lya polarization measurements for a giant Lya neb- 
ula at 2; ~ 2.656. We do not detect any significant linear polarization of the Lya 
emission: Piya = 2.6 ± 2.8% (corrected for statistical bias) within a single large 
aperture. The current data also do not show evidence for the radial polariza- 
tion gradient predicted by some theoretical models. These results rule out singly 
scattered Lya (e.g., from the nearby AGN) and may be inconsistent with some 
models of backscattering in a spherical outflow. However, the effects of seeing, 
diminished signal-to-noise ratio, and angle averaging within radial bins make it 
difficult to put strong constraints on the radial polarization profile. The current 
constraints may be consistent with higher density outfiow models, spherically 
symmetric infall models, photoionization by star formation within the nebula or 
the nearby AGN, resonant scattering, or non-spherically symmetric cold accre- 
tion (i.e., along filaments). Higher signal-to- noise ratio data probing to higher 
spatial resolution will allow us to harness the full diagnostic power of polarization 
observations in distinguishing between theoretical models of giant Lya nebulae. 

Subject headings: galaxies: formation — galaxies: high-redshift — galaxies: evo- 
lution — techniques: polarimetric 



^TABASCO Postdoctoral Fellow; Department of Physics, Broida Hall, Mail Code 9530, University of 
California, Santa Barbara, CA 93106; mkpresco@pliysics.ucsb.edu 

^Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721 

^National Science Foundation, 4201 Wilson Boulevard, Arlington, VA 22230, USA 

^National Optical Astronomy Observatory, 950 North Cherry Avenue, Tucson, AZ 85719 



-2- 



Introduction 



Very large (~100 kpc) gaseous nebulae, detectable in the high-redshift universe by their 
spatially extended, luminous Lya emission, have received substantial attention as the pos- 
sible sites of ongoing galaxy formation. The question of what powers these Lya nebulae 
(colloquially known as Lya 'blobs') has been a subject of keen interest since this class of 



objects was first discovered more than a decade ago (e.g.. iFrancis et al.lll996l : iFynbo et al. 



19991 : iKeel et al.lll999l : ISteidel et al.ll2000l ). because the origin of the Lya presumably holds 
valuable insight into the dominant physical process at work in these complex regions. Ob- 
servational studies have unc overed eviderice that some Lyg nebula e are powered by grav- 



itational cooling radiation ( Nilsson et al. 2006[ 
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Matsuda et al.l 120071 : iDey et al.l 120051 : iPrescott et al.l l2009l ) . Theoretical studies have ex- 
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tion (ITaniguchi fc Shiova 
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Taniguchi et al.l 1200 ll : lOiikstra fc Loebl l2009l : iGoerdt et al. 



Faucher-Giguere et al.l 120101 ). However, the uncertainty regarding the degree of ob- 



scuration and internal geometry coupled with large uncertainties in the modeling of Lya 
radiative transfer has left our understanding of the Lya nebula phenomenon confused and 
incomplete. Additional constraints are therefore needed to shed new light on the question 
of what powers the Lya emission. In this Letter, we present the first constraints on the Lya 
polarization in a giant Lya nebula. 

The observ able Lya polarization from an astrophysical cloud is governed by two impor- 
tant factors fsee lhee fc Ahnlll998l : Eroeb fc Rvbic"killl999l : [Rvbicki fc Loeblll999l : lDiikstra fc Loeb 
20081 ) . First, since Lya is a resonance line and the primary recombination line of the most 
abundant element in the universe, Lya radiative transfer is dominated in most contexts by 
resonant scattering. Lya photons are repeatedly absorbed and reemitted by neutral H atoms 
until they are either destroyed via absorption by dust grains or escape the cloud. Due to 
thermal motions of the absorbing atoms, photons execute a random walk in b oth frequency 
and physical space (e.g.. lHarringtonlll973l : lNeufeldlll990l : lLoeb fc Rybickilll999l ) and are most 
able to escape an optically thick cloud when they have scattered into the wing of the line 
profile. Kinematic offsets within the gas - e.g., due to the Hubble expansion, outflows, or 
infall - can help shift the Lya photons out of resonance and allow them to escape more 
readily. 

The second important consideration is the processes that might lead to linear polar- 
ization of the observed Lya emission. In pure Rayleigh scattering, the phase function is 
such that a photon is twice as likely to scatter forward or backward as it is to scatter at a 
90° angle, but the polarization is highest for photons that scatter at 90°. Therefore, light 
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that is initially unpolarized can be highly polarized by scattering off neutral atoms if scat- 
tered through large angles. In the case of Lya scattering, the phase function and degree 
of polarization qualitatively resemble the case of pure Rayleigh scattering, but the details 



must be calculated quan tum mechanically (e.g., iBrandt fc ChamberlainlllQSQt IStenfldll980 
Brasken fc Kyrolal 119981 ). A key result is that Lya scattering in the wing of the line profile 
yield s polarization that is three times higher than Lya scattering in the line core (jStenfio 
19801 ). Thus, those photons that are most likely to escape (those in the line wing) are also 
those with the highest polarization. If, after the photons have been polarized, however, they 
remain trapped in an optically thick cloud, they will resonantly scatter many more times. 
Since polarization is more likely to be destroyed than preserved in any subsequent scattering 
event, the highest polarization results from single scattering; as a general rule, the more 
photons scatter, the less polarization will be observed. Robust predictions for more complex 
scenarios require detailed simulations because the ionization state and optical thickness of 
the cloud, its kinematics, and its geometry will all affect the degree of polarization. 

A number of theoretical papers have modeled the expected level of polarization in the 
Lyo; lir ie under a variety of a strophysical scenarios. Focusing first on the pre-reionization 
epoch, iRybicki fc Loebl (119991 ) showed that a cosmological Lya-emitting source embedded 
within the neutral intergalactic medium (IGM) should produce high levels of polarization in 
the Lya emission line and that the polarization should increase with radius from the source. 
At small radii, where photons are scattered many times before they shift out of resonance and 
escape, the polarization is predicted to be low (~ 14%). Photons that reach large distances 
from the source, however, have typically been scattered in one large jump and are traveling 
roughly radially. The Hubble expansion has shifted them into the wing of the line and their 
last scattering event has put them on a course to the observer, i.e. ~90° scattering, resulting 
in higher polarization (~ 32 — 60%). 

In the post-reionization epoch, the situation is somewhat more complicated, and simu- 
lations have thus far focused on spherically symmetric geometries. iDijkstra fc Loebl (|2008[ ) 
showed that in the case of backscattered radiation in a spherical galactic outflow, the Lya 
polarization is predicted to increase strongly with radius to levels as high as ~ 40%, depend- 
ing on the assumed Hi column density. The high polarization occurs because the kinematic 
offset of the gas shifts the photons out of resonance and into the wing of the line, where 
the induced polarization is higher and where a substantial fraction are able to escape after 
only one scattering event. The radial increase stems from the fact that photons seen from 
larger radii have scattered through larger angles in order to reach the observer. In spheri- 
cally symmetric, collapsing, optically thick clouds, the situation can be similar to that of a 



spherica 



outflow, with polariza tion increasing as a function of radius to values as high as 



35% (IDijkstra fc Loebl l2008l ). While photons undergo multiple scatterings in this case. 
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they do so in the line wing where polarization is higher. The radial polarization profile 
is related to the fact that the Lya radiation field is increasingly anisotropic with radius; 
photons emerging from large radii were preferentially traveling radially outward prior to 
scattering, leading to larger scattering angles and higher polarization. Finally, in the case 
of resonant scattering in the IGM, the Lya photons scatter a greater number of times and 
do so in the core rathe r than in the wing of the line, yielding lower polarization (~ 2 — 7%; 



Dijkstra fc Loebll2008l ). Dust was assumed to be minimal in these calculations, but in the 
presence of significant amounts of dust the polarization signature could be higher because 
the more highly polarized photons are also the ones that scatter the least and are therefore 
more likely to escape before they are destroyed. Although these predicted polarization levels 
are encouraging from an observational perspective, it is important to remember that they 
were derived for the case of spherical symmetry and for specific choices of Hi column density 
and velocity profile; in more complicate d scenarios, e.g., cold a ccretion along filaments, the 



polarization may be substantially lower ( iDijkstra fc Loebll2009l ). The Lya polarization from 



these more complex scenarios has yet to be rigorously modeled. 

The existing theoretical predictions have thus far lacked an observational response. In 
this Letter, we present the first constraints on the Lya polarization within a giant Lya nebula. 
In Section [2l we discuss our observations and data reduction approach. Section |3] presents 
our results and discusses the implications of our finding that the measured polarization of 
the Lya emission in this Lya nebula is low (^5% overall or <12% at large radii). Our 
conclusions are given in Section |H 



Observations & Reductions 



We chose as our target a giant Lya nebula at z ^ 2.656 ( iDey et al.l l2005l . hereafter 



LABdOS). With a diameter of roughly 150 kpc and a Lya luminosity of ~ 10 erg s~ , it is 
one of the largest and most luminous Lya nebulae known, making it a prime candidate for 
follow-up observations. LABd05 was first discovered thanks to its strong Spitzer/MIPS 24/im 
emission, stemming from an associa ted obscured AG N that is offset by ~20 kpc (projected) 
from the peak of the Lya emission (IDey et al.ll2005l . Prescott et al. 2011, in preparation). 



We obtained imaging linear polarimetry centered on the Lya line using the B ok 2.3m 
Telescope and the CCD Imaging/Spectropolarimeter, SPOL (ISchmidt et al.lll992al ). during 
the nights of 2007 May 13-15 UT. The instrument is a dual-beam polarimeter built around 
a rotating semi-achromatic half-waveplate and a Wollaston prism. The SPOL field-of-view 
at the Bok Telescope is 100 x 100 pixels = 51" x 51", which for our observations contained 
both LABd05 and two nearby sources that allowed us to verify the position of each dithered 
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frame. 

Maximizing the sensitivity of these observations to polarization in the Lja line required 
a high throughput narrow-band filter centered at Lja at the redshift of the nebula. To 
achieve this, we selected a filter with an intrinsic central wavelength slightly to the red 
of Lya but with high throughput (A^ = 4481. 19A, FWHM = 56.18A, T^ax = 63.83%), 
and used the fact that tilting an interference filter relative to the incident beam has the 
effect of shifting the central wavelength to the blue. Using filter transmission measurements 
taken at two different angles, we computed the effective dielectric index of refraction of the 
filter {rii = 2.1) and used the standard relation (A = Ac[l — ^^^]^''^) to determine that an 
inclination angle of 6* = 15?4 would be needed to shift the central wavelength of the filter to 
match the Lya line at this redshift {XLya = 4446A). We designed an aluminum shim to tilt 
the filter in the SPOL filter holder by this angle during the observations. Tests done using 
the instrument as a filtered grating spectrometer at the beginning of the run confirmed that 
the filter bandpass was correctly centered on Lya aX z ^ 2.656. 

Over three nights we obtained 18 x 3200s = 16 hours of total integration and followed 
the standard procedure for polarimetric observations with dual-band polarimeters. The 
observations suffered from light to moderate cirrus and poor seeing of l'.'5-2'.'3, which in turn 
led to increased guiding errors, but the advantage of dual-beam polarimetric observations is 
that changes in transmission cancel during the reduction of the data. To derive the Q Stokes 
parameter, we obtained 2 separate Q integrations with initial A/2 waveplate positions of 0° 
and 45°, respectively, relative to the instrumental zeropoint. Each Q integration was divided 
into 4 individual exposures (200s each), one for each of 4 equivalent waveplate positions. 
After each individual 200s exposure, the shutter was closed, the waveplate was rotated by 
90°, and the shutter was reopened to continue the integration. To derive the U Stokes 
parameter, the same procedure was used to produce 2 separate U integrations but with the 
initial waveplate position set to 22?5 and 67?5, respectively. After each Q and U sequence 
(3200s), the telescope was dithered by ~ 1 — 2" in order to reduce the effects of flat-field 
errors on the total flux image. 

We reduced the data using IRAF. We subtracted the bias using the overscan region but 
did not apply a dark correction as the dark current was measured to be negligible (~0.3 
counts per pixel in a 300s exposure). Domefiats for each Stokes sequence were stacked 
and divided by the median value; the resulting fiatfield frame was applied to all images. 
Cosmic rays were removed from individual exposures using xzaplll The standard IRAF- 
based reduction package for SPOL {impaired, written by G. Schmidt) was applied to the 



1 Written by Mark Dickinson (NOAO). 
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image/waveplate sequences to yield the polarization results. 

We determined the instrumental polarization using observations of unpolarized standard 
stars (GD319, BD+33 2642, and BD+28 4211). As the instrumental polarization was neg- 
ligible (< 0.1%), we did not apply any correction to the data. To measure the polarimetric 
efficiency of the instrument, we observed a Nicol prism illuminated with the flatfield lamp 
through a pinhole aperture. This effectively simulated an observation of a star with 100% 
polarization and yielded a polarimetric efficiency measurement for the A/2 waveplate of 97% 
in this wavelength bandpass. The science data were corrected accordingly. Observations 
of interstellar polarization standard stars HD 155528 and Hiltner 960 were made to rotate 



the me asured polarization position angle to the equatorial coordinate system (ISchmidt et al. 



1992bl ). The dithered offsets between individual frames were calculated using a nearby source 



within the image and used to shift each exposure prior to generating the final image stacks. 



3. Results & Discussion 

The mean total flux image is shown along with the Q and U Stokes parameter images 
in Figure [H The Lya emission from LABd05 is clearly detected in the total flux image at 
a signal-to-noise ratio of 36 (16 pix=8'.'2 diameter aperture). In the case of polarization due 
to light scattered from the AGN (shown in Figure [T] as a red cross located at a position 
angle of 12?5 relative to the Lja peak), we would expect the electric vector position angle of 
polarization to be aligned roughly E-W across the nebula. We therefore start by measuring 
the Q and U Stokes parameters within a 16 pix (8'.'2) diameter circular aperture centered 
on the peak of the Lya emission and find that Pots = \/'Q^~+TP = 3.8±2.8%, with the 
uncertainty (cr) derived from photon statistics. As a check, we performed the same aperture 
photometry on the 18 individual (unstacked) Q and U frames and find that the quoted 
uncertainty yields indices of unity for the mean Stokes parameters derived from these 
18 independent measurements. When we correct for statistical bias - an important effect 
at low polarization signal-to-noise ratios - the polarization from the circular aperture is 



P ^ PnhJl — {(^Pobsl PobsTY^"^ = 2.6±2.8%, where we assume ap^^^ ^ a (j War die &: Kronberg 



m4 \. 



Our single aperture measurement of LABd05 is consistent with low Lya polarization 
(<5%). An idealized scenario in which Lya photons from the obscured AGN that lies to the 
north are scattering once off a nearby cloud is clearly not allowed by these observations. 

In the case of a single source of Lya photons centrally located within the cloud (e.g., 
as in the case of spherical outflow or infall), the polarization vectors would be oriented 
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perpendicular to the radial vector from the center of the nebula (i.e., in concentric rings) 
and would cancel when measured with a circular aperture covering the entire nebula. We 
therefore measured the Q and U Stokes parameters within subapertures (Figure [2]) and 
correct the individual subaperture estimates for the position angle of the radial vector to 
each subaperture. In Figure [2] we plot these rotated Q and U Stokes parameters along with 
the resulting polarization fraction averaged within bins in radius. Bins were chosen such 
that each bin has a signal-to-noise ratio of >5 in the total flux image, but we note that 
the innermost radial bin is compromised by poor seeing. We find that the polarization in 
each radial bin is consistent with zero and that there is no evidence for a radial polarization 
gradient. 

The lack of significant Lya polarization when measured in radial subapertures is difficult 
to interpret due to the effects of seeing, decreased signal-to-noise ratio, and angle averaging 
within individual bins. To ma ke a fair comparison, we created simulated observations based 



on the theoretical models of iDijkstra fc Loebl (120081 ): two models of backscattering in a 



spherical outflow with Hi column densities of 10^^ cm~^ and 10^° cm~^, and a model of 
a spherically symmetric infalling, optically thick cloud. The backscattered outflow model 
is described in terms of a parameter etc, corresponding to the outer edge of the outflow. 
For the purposes of comparison, we set ac = 5" (chosen such that the simulated peak Lya 
surface brightness matches that observed for LABdOS, after scaling the model appropriately 
by redshift and luminosity); however, our qualitative results are not sensitive to this choice. 
The simulated data were convolved with the typical seeing (FWHM = 2'.'0) and binned 
radially in the same way as the actual data. The simulated measurements are overplotted 
in Figure [21 

The observed upper limit on the polarization at the center of the last radial bin is 
< 11.5% (lo"). Given the effects of seeing and angle averaging, this corresponds to a Icr limit 
of ^25% on the intrinsic polarization at a radius of ~5". We flnd that the observed data are 
marginally inconsistent with the lower density (higher polarization) model of backscattered 
emission from a spherical outflow, but we cannot rule out higher density outflow scenarios or 
the case of a spherically symmetric infalling cloud or resonant scattering in the local IGM. 

There are two additional possibilities that are likely consistent with low Lja polariza- 
tion. The flrst possibility is Lya powered by photoionization from extended star formation 
or the AGN, which would be expected to either be unpolarized (in the optically thin case) 
or have low polarization similar to the resonant s catter ing scenario. The second is non- 



spherically symmetric cold flows: IDijkstra fc Loebl (120091 ) point out that Lya escaping from 



one fllament is not likely to be polarized by scattering off another due to the low predicted 
volume flUing factor of cold filaments. However, clear observational constraints on the cold 
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flow phenomenon are not yet available, and rigorous modeling of the expected Lya polariza- 
tion in cold flows has yet to be done. 

4. Summary 

In this Letter, we present the first observational constraints on the Lya polarization of 
spatially extended radio-quiet Lya nebulae. The polarization fraction of the Lya emission 
within LABdOS is 2.6±2.8% measured using a large circular aperture, which is inconsistent 
with a simple model of Lya photons (e.g., from the AGN) singly scattering off a nearby 
cloud. We see no evidence for a significant polarization gradient as a function of radius, 
a result that may rule out models predicting high polarization (~ 40%) from backscatter- 
ing in a spherical galactic outflow. However, the effects of seeing, decreased signal-to-noise, 
and angle averaging within radial bins make it difficult to constrain models that result in 
somewhat lower polarization: higher column density outflow models, spherically symmetric 
infall models, and models that invoke photoionization, resonant scattering, or cold flows with 
non-spherically symmetric geometries. These observations provide the flrst constraint on the 
Lya polarization of giant Lya nebulae, but as the radial polarization proflle is important 
for distinguishing between models, unleashing the full diagnostic power of polarization ob- 
servations of giant Lja nebulae will require higher signal-to-noise ratio data and analysis at 
higher spatial resolution. 
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Fig. 1. — Imaging polarimetry of LABdOS. The position of the Lja emission peak is indicated 
by the outer black crosshairs; the position of the obscured AGN is shown as a red cross. The 
lefthand panel shows the total flux image, and the middle and righthand panels give the 
Q and U Stokes parameter images. The spatially extended Lya emission is clearly visible 
in the total flux image (~ 36a in a 8'.'2 diameter aperture) but no polarization is detected 
(2.6±2.8%). 
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Fig. 2. — (Left) Subapertures used to measure the radial polarization profile are shown 
overlaid on the total flux image. The position of the obscured AGN is shown (red cross) 
at a position angle of 12?5 relative to the center of the Lya nebula. (Right) Rotated Q 
and U Stokes parameters and the resulting polarization fraction measured within radial bins 
using a weighted mean and corresponding error. Measurements for the innermost radial 
bin are compromised by poor seeing and are denoted with open circles. Simulated results 
based on theoretical models are overplotted (dotted lines). Models of backscattering in a 
spherical outflow - for column densities of Nhi = 10^^ cm~^ (red) and 10^" cm~^ (green) - 
and a spherical infall model (blue) are shown. Our observations put an upper limit on the 
observed polarization of < 12% (la) at large radii, corresponding to an intrinsic polarization 
limit of < 25%. 



